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ABSTRACT 

Helium-rich subdwarf B (He-sdB) stars represent a small group of low-mass hot stars with 
luminosities greater than those of conventional subdwarf B stars, and effective temperatures 
lower than those of subdwarf O stars. By measuring their surface chemistry, we aim to explore 
the connection between He-sdB stars, He-rich sdO stars and normal sdB stars. 

LS IV— 14°116 is a relatively intermediate He-sdB star, also known to be a photometric 
variable. High-resolution blue-optical spectroscopy was obtained with the Anglo- Australian 
Telescope. Analysis of the spectrum shows LSIV— 14°116 to have effective temperature 
T c ff = 34 000 ± 500K, surface gravity \ogg — 5.6 ± 0.2, and surface helium abundance 
"■He = 0.16 ± 0.03 by number. This places the star slightly above the standard extended hori- 
zontal branch, as represented by normal sdB stars. The magnesium and silicon abundances 
indicate the star to be metal poor relative to the Sun. 

A number of significant but unfamiliar absorption lines were identified as being due to 
germanium, strontium, yttrium and zirconium. After calculating oscillator strengths (for Ge, 
Y and Zr), the photospheric abundances of these elements were established to range from 
3 dex (Ge) to 4 dex (Sr, Y and Zr) above solar. The most likely explanation is that these 
overabundances are caused by radiatively-driven diffusion forming a chemical cloud layer in 
the photosphere. It is conjectured that this cloud formation could be mediated by a strong 
magnetic field. 
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1 INTRODUCTION 

Subdwarf B stars are low-mass core helium burning stars with ex- 
tremely thin hydrogen envelopes. They behave as helium main- 
sequence stars of roughly half a solar mass. Their atmospheres 
are generally helium deficient; radiative levitation and gravita- 
tional settling combine to make helium sink below the hydrogen- 
rich surface ( lHeberlll986h , to deplete other light elements, and 
to enhance abundances of heavy elements in the photosphere 
dO'Toole & Heberll2006l) . 

However, almost 5% of the total subdw arf population com- 
prise stars with helium -rich atmospheres dGreen et alj 1 19861 : 
I Ahmad & JefferyJ l2006h . Th ey have been va riously classified 
as sdOB, sdOC and sdOD dGreen et ail 1 19861) stars, bu t more 
recently as He-sdB and He-sdO stars dMoehler et ail Il99d : 
lAhmad & Jeffervll2004h . Their optical spectra are characterised by 
strong Hel (He-sdB) and Hell (He-sdO) lines. 



* E-mail: nas@arm.ac.uk 
f E-mail: csj@arm.ac.uk 



LSIV— 14° 1 16 was classified as an sdO star by 

iKilkennv & Pauls! dl990h and as an He-rich sdO star bv lViton et al.l 
1 1991 ). Jefferv et al.1 1 1996]) includ e d this star in their catalogue 
of He-sdB stars. lAhmad & Jeffervl (2005) reported pulsations in 
LS IV-14° 116 with periods of 1950 and 2900 s, the first discovery 
of pulsation in a He- sdB star. However, it is only mildly helium 
rich, with n Hc « 0.2 dAhmad & Jeffery||2003l) . 

This paper reports the first detailed abundance analysis of 
LSIV— 14° 116 using high resolution optical spectroscopy. In par- 
ticular, it reports the first discovery of strontium, germanium, yt- 
trium and zirconium in the optical spectrum of a hot subdwarf. 



2 OBSERVATIONS 

Three spectra of LSIV— 14° 116 were obtained with the Univer- 
sity College London Echelle Spectrograph (UCLES) on the Anglo- 
Australian Telescope (AAT) on 2005 August 26. Each exposure 
was 1800 s in duration. The observations w ere obtained, reduce d 
and analyzed in the same way as described bv lNaslimetai]d2010h . 
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Figure 1. YlII, ZrlV, SrII and Gem lines in LS IV- 14° 116. 

The AAT/UCLES spectrum of LSIV-14°116 shows Nil, 
NlII, ClI, Cm, Oil, Sim, and Mgll lines, together with strong Hel 
and HeIlA4686. Unlike other He-sdB stars, the spectrum displays 
relatively strong hydrogen Balmer lines. While analysing this spec- 
trum we noted a number of strong unidentified lines. Using the 
NIST database of atomic spectra, we discovered that these lines 
are due to ZrlV, YlII, Gem and SrII (Fig. [T). As far as we know, 
ZrlV, YlII, and Gem have not been identified in any other ground- 
based astronomical spectrum, although Zr and Ge have been ob- 
serve d with large overabundances in ultraviolet spectra of sdB 
stars dO'Tooldl2004)O'Toole & Hebenl2006l : IChaver et alj|2006t 
iBlanchette et alJ l2008 ) . The SrII resonance lines are seen in chem- 
ically peculiar late-B stars (HgMn stars, for example), but are 
completely unknown at late-O or early-B spectral types. These 
line identifications immediately pointed to a remarkable chemical 
anomaly in this star. However, the anomaly could not be measured 
immediately, since no oscillator strengths existed for the optical 
ZrlV, YlII, and Gem lines. 

We have not identified lines for any iron-group elements in 
this spectrum. There remain a small number of lines which we 
have been unable to identify. The most notable have wavelengths 
4007.40A and 4216.20A with equivalent widths 40 and 34 mA, 
respectively. 
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Figure 2. Theoretical line profile fits to Hel 4471 and H 7 in LS IV- 14° 1 16 
for model atmospheres with njj c = 0.100, T e g = 34 000 K and log g = 
4.5(0.5)6.0. 
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Figure 3. The loci of ionization equilibria for Cll/ClII, Nll/NlII and Hell, 
the profile fits to Hel and H lines, and the adopted solution. 



3 PHYSICAL PARAMETERS OF LS IV 14° 116 

We measured the effective temperature T e s, surface gravity log g, 
and elemental abunda nces of LS IV— 14° 116 using the same 
method as described bv lNaslim et alj feOlOh . We found the micro- 
turbulent velocity v t = 10 km s _1 by minimising the dispersion on 
carbon abundance measured from six Cm absorption lines; adopt- 
ing Vt = 5kms~ increases this dispersion by 25%. The grids of 
model atmospheres adopted for analysis were computed with 1/10 
solar metallicity, relative helium abundances rm c ~ 0.10, 0.20 and 
0.50 by number, and microturbulent velocity v t = 10 km s _1 . 

T c g and log g were measured using ionisation equilibrium and 
line profile fitting. Using appropriate model atmospheres, the ioni- 
sation equilibrium was established by balancing the abundances de- 
termined from Nil and NlII and from ClI and Cm lines in the spec- 
trum, as well as by fitting the equivalent width of the temperature 



sensitive Hell 4686 line. The surface gravity was established by fit- 
ting theoretical profiles to H/3,7, and 5, and to Hel 4471 (Fig. [2j. 
The coincidence of profile fits and ionisation equilibria was used to 
determine the overall solution illustrated in Fig. [3] 

We also measured the physical parameters T e ff, log<? and 
helium abundance rinc using the package SFIT which finds the 
best-fit solution within a grid of synthetic spectra. The model 
grid was defined with T cft = 32 000 (2000) 38 000 K, log g = 
4.50(0.5)6.00, n Hc = 0.10, 0.20, 0.30, 0.50 and v t = lOkms" 1 . 



Table 1. Atmospheric parameters 



r eff (K) 


log 9 


"He 


Source 


34 000 ± 500 


5.6 ±0.1 


0.160 ±0.03 


SFIT 


33 860 ± 1 800 


5.54 ±0.2 




Ion eq 


33 000 ±1000 


5.8 ±0.2 


0.2 


Viton et al. (1991) 


32 500 ±150 


5.4 ±0.1 


0.21 


Ahmad & Jeffery (2003) 
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Table 2. Individual line equivalent widths W\ and abundances ei 
adopted oscillator strengths gf, forLSIV-14°116. 
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Figure 4. The location of LS IV — 14° 116 on the log , q - T eft diagram, 
compared with normal sdB stars I Edelmann et al 1 l200i helium-rich sdB 
stars lAhmad & Jeffery||2003l : iNaslim et al]|2oiol) helium-rich sdO stars 
IStroeer et alj|2007l) . and JL87 lAhmad et al.ll2007l) . A zero-age main se- 
quence (Y = 0.28, Z = 0.02), a zero-age horizontal branch (M c = 
O.485M , Y c = 0.28, Z = 0.02) and a helium main sequence (Z = 
0.02) are also shown. 



The observed spectrum of LS IV— 14° 116 is shown along with the 
best fit model spectrum in Fig. [5] 

The atmospheric properties of LS IV— 14° 116 deduced from 
these measurements are given in Table Q] The projected rota- 
tional line broadening measured from Znv lines 4198.26A and 
4317.08A is vsini < 2kms . SFIT and ionization equilib- 
rium results agree to within the formal errors. Although all ioni- 
sation equilibria and profile fits do not converge on a single point, 
the intersections lie within 3200 K in T e a and 0.5 in logg. We 
took a weighted mean of these intersections to determine the ion- 
ization T c ff and logy. For subsequent abundance measurements, 
we adopted the grid model atmosphere with T e g — 34 000 K, 
loggr = 5.5, nH C = 0.20 and v t = lOkms -1 . 

In view of the atmospheric abundances reported below, the 
above procedure was repeated with a grid of models in which the 
abundances of all metals with atomic number ^ 20 (calcium) ex- 
cept iron were increased to 10 times the solar value. This has proved 
to be a successful strategy for the analysis of normal sdB star spec- 
tra JO'Toole & HebeiT2 006; Pereira 2010). However in this case, a) 
we failed to find any agreement between HeII4686 and the Nll/lII 
ionisation equilibrium, b) we found the depth of the theoretical H 
and Hel lines to be consistently deeper than observed and c) whilst 
leading to a 2000 K reduction in T e g , the new models also led to a 
significant increase in the measured strontium abundance. Without 
further diagnostics of high-Z abundances and chemical stratifica- 
tion, the 1/10 solar metallicity models were deemed satisfactory 
for the present paper. 



4 ATOMIC DATA 

In order to measure elemental abundances for all visible species, 
we measured equivalent widths of lines of Cll.Clll, Nil, NlII, Oil, 
Mgll, Sim, SrII, Gem, YlII and Zrm. Atomic data and oscillator 
strengths for most of these lines were available in our own compi- 
lation (see Tabled, but data for Sr, Ge, Y and Zr were not. 

The atomic data for the relevant transitions of these ions are 
displayed in Table [3] The labels of the upper and lower states of 



Ion 








A/A 


logs/ 


W\/m.} 




ClI 








4074.48 


0.204 


1 42 


8.20 


4074.52 


0.408 






4074.85 


0.593 


j 




4075.85 


0.756 


57 


8.08 


4075.94 


-0.076 






4267.02 


0.559 


130 


7.80 


4267.27 


0.734 






Cm 








4067.94 


0.827 


60 


7.70 


4068.91 


0.945 


59 


7.69 


4070.26 


1.037 


82 


8.37 


4186.90 


0.924 


76 


8.34 


4647.42 


0.072 


130 


8.08 


4650.25 


-0.149 


98 


8.01 


4651.47 


-0.625 


54 


7.99 


Nil 








3995.00 


0.225 


47 


7.65 


4041.31 


0.830 


45 


7.82 


4043.53 


0.714 


46 


7.95 


4236.86 


0.396 


1 73 


8.19 


4236.98 


0.567 


J 




4241.78 


0.728 


] 35 


7.85 


4241.79 


0.710 


J 




4447.03 


0.238 


48 


8.04 


4601.48 


-0.385 


35 


8.15 


4607.16 


-0.483 


42 


8.35 


4630.54 


0.093 


50 


7.88 


Nm 








4640.64 


0.140 


77 


8.19 


4634.14 


-0.108 


51 


8.13 


Oil 








4072.15 


0.545 


21 


7.42 


4414.90 


0.210 


33 


7.63 


4416.97 


-0.041 


26 


7.75 


Mgll 








4481.13 


0.568 


25 


6.85 


Sim 








4552.62 


0.283 


29 


6.32 


gf values: ClI 


Yan et al. 


11987), Cm 


Hibbert 


Hardorp & Scholz ( 1970); 


Bockasten 




Becker & Butler 


(1989), 


NlII iButlei 


(1984) 


Becker & Butler 


(1988). Sim iBecker & ButleJ 



(1966 



Nil 
On 



each transition include, in addition to the angular momenta and 
parity, the electron occupancy of the valence subshells, as given 
by the Hartree-Fock approximation. In practice, the Hartree-Fock 
(HF) approximation is not sufficiently adequate for th e calculation 
of osc illator strengths, as demonstrated by the work of iBrage et al.l 
dl998h for SrII. Instead, each wave function needs to be represented 
by a configuration interaction (CI) expansion, though the dominant 
(largest CI mixing coefficient) remains the HF configuration in all 
the cases considered here, and therefore constitutes an appropriate 
part of the label. 

For the SrII tr ansitions, w e have merely quoted the oscillator 
strength results of IBrage et alj ([1998), since their work is a sub- 
stantial MCHF calculation and their final r esults agree well with 
the experimental studies of iGrevesse e t al. ( 199 % and even with 
the conceptually simpler calculation of Migdalek & B avliss! j 19791) 
who represented core polarisation by model potentials rather than 
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Table 3. Atomic data for optical lines of Sr, Ge, Y and Zr, and abundances derived for LS IV— 14°116. 



A/A 


Confi 


guration 


Ei/eV 


Ref 


log gf 


Ref 






SrII 


















4077.71 


5s 2 S 1/2 


- 5p 2 P 3 /2 


0.00 


1 


0.142 


1 


24 


7.00 


4215.52 


5s 2 S 1/2 


- 5p 2 P 1/2 


0.00 




-0.175 




10 


6.91 


Gem 


















4178.96 


4s5s 3 Si 


— 4s5p 3 P 2 


19.66 


2 


0.341 


5 


63 


6.36 


4260.85 


4s5s 3 Si 


— 4s5p 3 Pi 


19.66 




0.108 




40 


6.34 


4291.71 


4s5s 3 Si 


— 4s5p 3 Po 


19.66 




-0.368 




10 


6.14 


YlII 


















4039.602 


4/ 2 P 7 /2 


- 5g 2 G g/2 


12.53 


3 


1.005 ] 


5 


21 


6.12 


4039.602 


4/ 2 P 7 /2 


- 5g 2 G r / 2 


12.53 




-0.538 J 








4040.112 


4/ 2 P 5/2 


- 5g 2 G T/2 


12.53 




0.892 




21 


6.24 



ZrlV 

4137.435 5d 2 D 3/2 - 6p 2 P 3 ° /2 18.18 4 -0.625 5 10 6.20 

4198.265 5d 2 D 5/2 - 6p 2 P 3 ° /2 18.23 0.323 86 6.43 

4317.081 5d 2 D 3/2 - 6p 2 P° /2 18.18 0.069 64 6.49 

4569.247 5g 2 G 7/2 - 6h 2 H° /2 25.65 1.127 ] 101 6.76 

4569.247 5g 2 G 9/2 - 6h 2 H° 1/2 25.65 1.216 J 

References :!. Brage et al. 11998). 2. Lang~1929). 3. Epstein & Reader Il975l) . 4. 

Reader & Acquista Ul997h , 5. this paper. 



by explicit CI. However, transition data for the other ions has not 
been calculated to such accuracy, and therefore we undertook to 
provide that data ourselves. 

For ions with valence electrons outside a number of closed 
subshells, there are three types of electron correlation which im- 
prove upon the HF approximation: valence shell correlation (with 
the core described as in the HF approximation); core- valence cor- 
relation (or core polarisation), in which the influence of the core on 
the valence electrons is included; core-core correlation. For mono- 
valent ions, only the second and third of these effects occurs, in 
accordance with the work of iBrage et al] il998h . We describe be- 
low our work on the three ions Gem, YlII, ZrlV. 



vice versa. Core-core correlation is likely to have a very small ef- 
fect on the oscillator strengths. 

Our oscillator strengths (length form) are shown in Table [3] 
and are obtained using experimental wavelengths. The calculated 
energies agree with experiment to within 2%. The length and ve- 
locity forms also agree to within 5%. We also undertook the cor- 
responding calculations in LS coupling and found approximately 
the same level of accuracy. We noted that, unlike the SrII work, 
core polarisation had only a small effect. This is partly because 
the orbitals of the transition (5s and bp) are further from the core, 
and partly because the effective core experienced by these electrons 
also includes the 4s electron. 



4.1 Ge in 

This is the only ion considered here which has two valence elec- 
trons and therefore the only ion in which valence shell correlation 
must be included. Much work has bee n published on tran sitions 
amongst the n = 4 levels (for example. IChen & Chendl201ol) . but 
almost nothing on transitions involving higher levels. We under- 
took our calculations o f oscillator strengths using the confi guration 
interaction code CIV3 faibbertl 19751 ; iHibbert et al J 199 lh . 

Relativistic effects, giving rise to fine-structure splitting, were 
included using the Breit-Pauli approximation, but all orbital opti- 
misations were carried out in LS coupling. The radial functions of 
the orbitals are expressed in analytic form as sums of Slater-type 
orbitals. We used (in the absence of functions f or Gem) the radial 
functio ns of the ground state of Gell given by lOementi & RoettH 
(1974) but reoptimised the outer orbital functions 4s and 4p on the 
energy of the 4s4p 3 P° state of Gem. With these orbitals fixed, 
we optimised 5s and 5p on the energies of 4s5s 3 S and 4s5p 3 P°. 
Subsequent to this, we optimised 4d and 4/ on the energy of 4s5p 
3 P° and 5d on 4s5s 3 S. Valence-shell correlation was represented 
by constructing all possible configurations using these orbitals, but 
keeping a common core of filled subshells up to 3d. 

Core-valence correlation was represented by configurations 
with a single electron removed from 3d 10 and replaced by Qp or 
5/, optimised on 4s5s 3 S, with 4s5s replaced by 4s4p/4s5p and 



4.2 Y III 

We proceeded in a similar manner for YlII: the starting point were 
the radial functions for the ground state of Yll given by Clementi 
and Roetti (1974), and then reoptimised 4s, 4p along with 4/ on 
the 4s 2 4p 6 4f 2 F° state, and then the 5g orbital function on the 
2 G state. In this case, no valence-shell correlation is possible, so 
core polarisation is the dominant correlation effect modelled ex- 
plicitly by opening the 4p and 4s subshells. Specifically, we added 
to 4s 2 4p 6 5g 2 G the configurations 4s 2 4p 5 (4d6ft + 4d4/ + 4d5/) 
and 4s4p 6 (4/5p + 5p5/ + 5p6h), with similar additions for the 
2 F° state. We also added limited core-core correlation by means of 
configurations in which 4s is replaced by 4d in the main configura- 
tions, as well as 4p 6 by 4p 4 4d 2 . 

The resulting oscillator strengths are shown in Table[3] Again, 
calculated transition energies agree to with 4% of experiment. We 
consider that the oscillator strengths are accurate to better than 5%. 
As with Gem, the effect of core polarisation is quite small, which 
is not surprising in view of the higher orbital angular momenta of 
the outer electron. 

4.3 Zr IV 

A similar process to that adopted for YlII was carried out for ZrlV, 
with the core orbitals chosen as those of the Zrll ground state, while 
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Table 4. Mean abundances & for LS IV- 14° 116. 



5 ABUNDANCES 



El. 


LSIV-14°116 


sdB 1 - 5 


JL 87 b 


He-sdB Y 


Sun 8 


H 


11.83 ± 0.07 


12.0 


11.6 


< 8.5 - 11.1 


12.00 


He 


11.15 ± 0.05 


7.9-11.4 


11.3 


11.5 


10.93 


C 


8.04 ±0.22 


6.9-7.5 


8.8 


7.1-9.0 


8.52 


N 


8.02 ± 0.20 


7.4-7.6 


8.8 


8.3-8.7 


7.92 


O 


7.60 ±0.17 


6.5-8.5 


8.6 


7.1-7.3 


8.83 


Ne 


< 7.6 


6.5-8.5 


8.31 


8.2-8.5 


[8.08] 


Mg 


6.85 ± 0.10 


5.6-7.6 


7.4 


7.1-8.3 


7.58 


Si 


6.32 ± 0.12 


5.5-7.7 


7.2 


6.8-7.4 


7.55 


Ar 


< 6.5 


6.0-9.0 






[6.40] 


Sc 


< 5.3 


5.0-7.0 






3.17 


Ti 


< 6.0 


5.3-9.0 






5.02 


V 


< 6.5 


6.0-8.5 






4.00 


Cr 


< 7.0 


5.5-8.0 






5.67 


Fe 


< 6.8 


6.5-8.1 






7.50 


Ge 


6.28 ±0.12 


1.5-5.5 






3.41 


Sr 


6.96 ±0.15 








2.97 


Y 


6.16 ±0.10 








2.24 


Zr 


6.53 ±0.24 


2.7-^4.9 






2.60 



l. lEde lmannetal. 12003), 2. Geier et al. 1 200 8). 3.IGeier et all <2 010). 
4. lO'Toole & HebeJ fe006l). 5. IChaver et alj 42006h 6. lAhmad et all 
( 2007), 7. Naslim et all told). 8,lGrevesse & Sauvall jl99^) 



Table 5. Abundance errors <5e; due to representative errors in T a g, logg 
and vt 



Element 


8T cff = 1000 K 


<51ogg = 0.2 


8v t = 5 


C 


+0.05 


+0.03 


-0.11 


N 


-0.03 


+0.02 


-0.08 


O 


-0.08 


-0.004 


-0.04 


Mg 


-0.07 


-0.01 


-0.03 


Si 


-0.13 


-0.02 


-0.04 


Ge 


-0.13 


-0.02 


-0.05 


Si 


-0.24 


-0.08 


-0.01 


Y 


-0.08 


-0.02 


-0.03 


Zr 


+0.02 


+0.02 


-0.11 



4s, Ap were reoptimised along with 5g on the ZrlV configuration 
5g 2 G and, with these core orbitals, 6h was optimised on the 6h 
2 H° state. In a similar manner to Yin, orbitals 4d, 5s, 5p were 
introduced to allow for core polarisation. Two other orbitals, 4/ 
and 7i, were introduced to allow for a more angular flexibility in 
the outer electrons. Again the effects of core polarisation are small, 
and length and velocity forms of the oscillator strengths agreed to 
around 3%. 

A separate calculation was undertaken for the 5d-6p tran- 
sitions. We reoptimised 4p along with Ad on the 2 D state, and 
then with these functions optimised 5d, 5p and 6p on their respec- 
tive states. To represent core- valence correlation, additional orbitals 
5s, 7p, 6d, 4/ were optimised to allow for configurations with one 
electron removed from the 4p subshell and 8p in order to polarise 
the 4s subshell. The calculated transition energies agreed very well 
with experiment (to within 1%). We estimate that the oscillator 
strengths are correct to within better than 10%. The effect of core 
polarisation was still small, principally due to the high n-values of 
the outer electrons. 



Individual line equivalent widths and abundances are given in Ta- 
bles [3] and flfl Abundances are given in the form e, = log m + c 
where log Sia^n,: = logEiaiWi© = 12.15 and m are atomic 
weights. This form conserves values of for elements whose abun- 
dances do not change, even when the mean atomic mass of the mix- 
ture changes substantially. 

Mean abundances for each element are given in Table [4] The 
errors given in Table [4] are based on the standard deviation of the 
line abundances about the mean or the estimated error in the equiv- 
alent width measurement. Systematic shifts attributable to errors in 
T e ff, log g and v t are given in Table [5] The final best-fit spectrum 
using the adopted best-fit model and with the elemental abundances 
from Table|4]is shown in Fig. [5] together with identifications for all 
of the absorption lines in the model. 

Table|4]also shows a representative range of abundances mea- 
sured for normal sdB stars, for the intermediate helium sdB star 
JL 87, for helium-rich sdB stars and the Sun. In common with many 
helium-rich sdB stars, LS IV— 14°116 appears to be mildly metal- 
poor (Mg and Si are assumed as proxies for the mean metallicity). 
Relative to the Sun, nitrogen is unremarakble, but relative to Mg 
and Si, LS IV- 14° 116 is arguably N-rich. 

With a helium-to-hydrogen ratio of 0.2, LS IV-14°116 falls 
well between the normal sdB stars, in which heliu m is strongly 
deple t ed, and the truly hel ium-rich sdB and sdO stars dNaslim et al.1 
l20irAIStroeeretal.ll2007l). Only JL 87, with n H /n He « 1 might 
be comparable ( lAhmad et al]2007l) . However, JL 87 is strongly C- 
and N-rich, and shows no evidence for exotic species such as those 
described here. 

Heavier elements, including argon and the iron-group ele- 
ments scandium, titanium, vanadium, and chromium are frequently 
seen to be overa bundant in H-ri ch sdB stars with higher effective 
temperatures dGeier et al.1 120101) . LS IV— 14° 116 does not show 
any detectable lines due to these elements. We estimated upper limit 
abundances assuming a minimum observable equivalent width of 
5mA (Table |4j. Overabundances of Sc, Ti, V and Cr cannot be 
ruled out. 

What is undeniable is that LS IV— 14° 116 shows overabun- 
dances of Ge, Sr, Y and Zr which are unprecedented in any hot 
subdwarf and, possibly, in any other star (Sr, Y and Zr are only ~ 
2 dex overabundant in Przybylski's star HP 1 01065, although heav - 
ier elements do show 4 dex overabundances: ICowlev et all faoO); 
IShulvak et aflfeOlOh ). 



6 DISCUSSION 

Both normal sdB stars and helium-rich subdwarfs are chemically 
peculiar. The former are depleted in helium and in other ele- 
ments wi th Z < 20, but enriched in all elements Z > 20 ex- 
cept iron dO'Toole & Heberll2006| ; |Geier et alj|20ld) . The consen- 
sus explanation is that chemical diffusion operates in an atmo- 
sphere where gravitational settling and radiative levitation com- 
pete to elevate ions to layers in which their specific opacities are 
maximised. Abundanc e informatio n for helium-rich subdwarfs is 
less extensive dNaslim et al . 2010), but the absence of hydrogen 
and the enrichment of nitrogen and (sometimes) carbon point to 
nuclear-processed helium dominating the surface composition. A 



1 Abundances were calculated with i mproved model atmo spheres, which 
give values different to those shown bv lNaslim et aljfeoiOh 
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Figure 5. The merged AAT/UCLES spectrum of LSIV-14°116 along with the best FIT model of T eff = 34 000 K and log g = 5.50. Abundances areas in 
Tableg] 



few intermediate-helium sdB stars (with 0.1 < n-Hc < 0.9) may 
be transition objects; if helium-rich sdBs are contracting towards 
the extended horizontal branch, intermediate-helium sdB stars may 
represent those stars in which diffusion has been effective suffi- 
ciently long to depl ete the surface helium by a significant amount 
dNaslim et al.ll2010h . 

6.1 Diffusion 

However, chemical peculiarities of the magnitude seen on the sur- 
face of LS IV— 14° 116 demand a more sophisticated explanation. 



Are they produced by some diffusion process which has led to ex- 
treme concentrations (i.e. clouds) in the line-forming region of the 
photosphere? Or are they a consequence of dredging-up heavily nu- 
clear processed material from the intershell region of an AGB star? 

LSIV— 14°116 would barely qualify as "helium enriched" 
were it a main-sequence star. Only because sdB stars are generally 
so helium poor was it considered unusual in the first place. So the 
question of whether diffusion is responsible for other abundance 
anomalies is legitimate. 

The elements Sr, Y, and Zr are recognized as products of s- 
process neutron-capture reactions which occur in the intershell in 
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asymptotic giant-branch (AGB) stars. They can appear on stel- 
lar surfaces following either third dredge- up on the AGB, a very 
late th er mal pulse (e.g. FGSge, V4334 Seei ljefferv & SchonberneJ 
(2006); Asplun d et al] dl999l) ). or as a co nsequence of a whi te 
dwarf merger, (e.g. V1920Cyg, HD124448: [Pandev et al.l d2004h ). 
The difficulties here are: (a) that none of these processes has pre- 
viously been indicated to produce such a large excess; 1-2 dex is 
more common, (b) the absence of a characteristic carbon and/or 
oxygen excess, and (c) there is no evidence that LS IV— 14° 116 is 
a post-AGB star in any sense. 

Some abundances anomalies have been recognized as due to 
pollution from an external source. Dust-fractionated accretion from 
the interstellar medium tends to produce modest anomalies which 
are correlated with condensation temperature. This is not the case 
here. Accretion from the wind of a now unseen AGB companion 
can also be ruled out by the absence of a carbon and/or oxygen 
excess. 

The principal advantage of the diffusion hypothesis is that it 
allows for a very high overabundance of a given species in a very 
narrow layer of the stellar atmosphere, but without requiring an ad- 
ditional source of that material providing that there is a commensu- 
rate depletion in other layers. Since the line-forming region of the 
stellar atmosphere (optical depth 0.001 < r < 0.1) contains ~ 4± 
lxl0~ 14 M o compared with an estimated 7 x 10" 13 M Q for the to- 
tal mass of the atmosphere (r < 10) and ~ O.OOIM0 for the mass 
of the stellar envelope, the observed excess can be easily under- 
stood. It is important to note that overabundances of strontium and 
yttrium of up to 3 dex are well-documented in c hemically-peculiar 
A (Ap) and mercury-manganese (HgMn) stars {Cowley & Hubrid 
120081; ICowlev et alfeoid : iDworetskv et alj|2008l) . where they are 
commonly attributed to radiatively-driven diffusion producing a 
strongly stratified atmosphere. 



sitive to probe the behaviour of Ge, Sr, Y, and Zr lines, as well as 
the H and He lines would be indicative of differential motion within 
the stellar atmosphere. 

The second possible link refers to observations of intermediate 
helium stars on the main sequence, so metimes known as B p(He) 
stars. The class prototype is a OriE dGreenstein & WallersteirJ 
1958), a B2Vp star with a 1.1 d variation in light, a 
chemically-inhomogeneous s urface and a strong magnetic field 
dPedersen & Thomsenlll977l ; iThomsenl fl974l ; lLandstreet & Borral 
1 19780 - Strong magnetic fields are frequentl y associated with large 
abundance anomalies i n early-type stars (Bohlender et al.l Il987l ; 
lHunger&Grootell999h . In simple terms, a magnetic field radically 
alters the opacity of a given ion. Specifically, by splitting energy- 
level degeneracies, Zeeman splitting increases the line opacity in 
a direction parallel to the magnetic field lines. Consequently, ra- 
diative levitation preferentially leads to the greatest accummula- 
tion of susceptible elements in regions where the magnetic field is 
orthogonal to the stellar surface, i.e. at the magnetic poles. Such 
anoma lies can be further exa cerbated by fractionation in the stellar 
wind dHunger & Grootdl 19991) . Strong che mical anomalies appe ar 
as "spots" and can lead to light variability dTownsend et al.l r2005). 

We are not currently suggesting that LS IV— 14° 116 is aclass- 
cial Bp(He) star. Its surface gravity is too high and so a connec- 
tion with the low-mass sdB stars seems more probable. However, 
we are suggesting that the same magnetically-driven physics could 
be responsible for the chemical anomalies. This suggestion can be 
tested. It requires that LS IV— 14° 116 have a strong magnetic field, 
and it is possible that the light variations could be associated with 
chemical variability. High-resolution time-resolved spectroscopy 
would indicate whether the surface is chemically homogeneous. 
Simultaneous photometry would demostrate any link to light vari- 
ability. Spectropolarimetry would indicate a magnetic field. 



6.2 Variability 

LS IV— 14° 116 i s know n to be a low-amplitude light variable. 
lAhmad & Jeffervl d2005l) report two periods (1950 and 2900 s) and 
suggest these are due to non-radial (/-mode oscillations of high ra- 
dial order, since the periods are too long to be p-modes. A dif- 
ficulty with this proposal is that q-modes are not predicted i n 
sdB stars hotter than about 29 000 K djefferv & SchonberneJioO^) : 
with r e a = 34 000 K, L S IV- 14° 116 sits in t he middle of the sdB 
p-mode instability zone ( Char pinet et al.l200ll) . Our three consecu- 
tive AAT spectra show no evidence of variability in line equivalent 
widths or radial velocity. This is probably a consequence of the low 
amplitude of the variability and the exposure times (1800 s) being 
long compared with the photometric periods. 

There are two possible links between photometric variabil- 
ity and unusual chemical composition. The first regards the blue 
edge of the g-mode instability strip, which is very sensitive to the 
metal abundance in the driving zone dJefferv & Saicll2006al) . Re- 
ducing the hydrogen a bundance also tends to destabilize pulsation 
(Jeffer v~& Saiol2006bl) but, in this case, the hydrogen abundance is 
not sufficiently depressed to make a major difference. Should the 
overabundances seen in Ge, Sr, Y and Zr be reflected in other el- 
ements and also be continuous from the atmosphere to the driving 
zone at ~ 2 x 10 5 K then, in all probability, there would be con- 
sequences for pulsational stability. However, as discussed above, 
there is no obvious source for such an excess throughout the stellar 
envelope. A strong requirement to demonstrate that the light vari- 
ability of LS IV— 14° 116 is due to pulsation is to investigate its 
radial-velocity behaviour. Observations which are sufficiently sen- 



7 CONCLUSION 

LS IV— 14° 116 is an intermediate helium-rich subdwarf B star, 
with a gravity slightly lower than that of normal sdB stars, and a 
surface helium abundance of about 16% by number. Taking sili- 
con and magnesium as proxies for the mean metal abundance, it 
is slightly metal poor (-0.8 dex) relative to the Sun. Since nitrogen 
is approximately solar, and carbon and oxygen are 0.4 and 1 . 1 dex 
subsolar respectively, the excess helium is overabundant, it is pos- 
sible that the surface consists of a mixture of hydrogen with CNO- 
processed helium. How diffusion affects C, N and O is not clear. 
Neither is it currently known whether LSIV— 14°116 is a single 
star or has an undetected binary companion. 

I n conjunction w ith observations of other helium-rich sdB 
stars dNaslim et al]|2010l) . it might be argued that LSIV- 14° 116 
was formed from either a helium white dwarf merger, or a late 
helium-flash episode, in which such a mixture could arise. Whether 
the observed helium/hydrogen ratio represents the relative contri- 
butions of the original material is difficult to establish; radiative 
levitation may already have caused a substantial fraction of the he- 
lium to sink as the star contracts towards the zero-age extended 
horizontal branch. 

What makes LSIV— 14° 116 remarkable, and distinct from 
any other hot subdwarf, whether helium-rich or not, is the pres- 
ence of lines due to germanium (Gem), strontium (Sill), yttrium 
(Yin) and zirconium (Zriv) in its optical spectrum. Measurements 
of these lines translate into overabundances of up to four orders 
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of magnitude of these elements in the line-forming region of the 
photosphere. 

Since a nuclear origin for such overabundances seems un- 
likely, it is argued that they arise due to radiatively-driven diffu- 
sion, where elements accumulate in layers of high specific opacity 
and where radiative and gravitational forces are in equilibrium. It 
is conjectured that a strong magnetic field might be responsible for 
the extreme overabundances observed in this case. How the unusual 
chemistry is linked to the photometric variability remains to be es- 
tablished. 
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